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Soliton pinning by long-range order in aperiodic systems 
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We investigate propagation of a kink soliton along inhomogeneous chains with two different con- 
stituents, arranged either periodically, aperiodically, or randomly. For the discrete sine-Gordon 
equation and the Fibonacci and Thue-Morse chains taken as examples, we have found that the 
phenomenology of aperiodic systems is very peculiar: On the one hand, they exhibit soliton pinning 
as in the random chain, although the depinning forces are clearly smaller. In addition, solitons are 
seen to propagate differently in the aperiodic chains than on periodic chains with large unit cells, 
given by approximations to the full aperiodic sequence. We show that most of these phenomena 
can be understood by means of simple collective coordinate arguments, with the exception of long 
range order effects. In the conclusion we comment on the interesting implications that our work 
could bring about in the field of solitons in molecular (e.g., DNA) chains. 
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The subtle interplay between nonlinearity and disor- 
der is being laboriously unveiled throughout the past few 
years |Q . A rich diversity of phenomena stems from such 
interaction, their manifestations being found in a num- 
ber of systems ranging from condensed matter physics 
to biophysics [g. A number of models have been set 
forth which capture the essential ingredients of those sys- 
tems while enjoying a canonical, non-specific view of the 
problem. Among the most successful of these models, 
the sine-Gordon (SG) equation is particularly remark- 
able both for its range of applicability and the possibil- 
ities it opens for study either in continuous or discrete 
version. Some of the physical situations well modeled by 
this equation are, for instance, Josephson junctions jl, 



Josephson junction arrays (J J A) |t||J, or DNA promoter 
dynamics Importantly, many realistic systems like 

DNA chains are neither periodic nor random, being in- 
herently close to quasi-periodic or aperiodic systems, so 
that the effects of long-range order may change the dy- 
namics of nonlinear excitations. 

In this Rapid Communication we concern ourselves 
with the problem of the behavior of kink solitons on 
lattices consisting of two different components, thereby 
focusing on issues inherently discrete similar to those 



of DNA or JJA dynamics. Our main aim here is to 
learn about the phenomenology of soliton propagation as 
a function of the order of the underlying lattice. We con- 
sider three main possibilities for our binary chain: pe- 
riodic, aperiodic, and random, which represent, respec- 
tively, full order, long-range order, and pure disorder. 
We show in the following that, while the periodic lattice 
exhibits basically the same features as the homogeneous 
case, the two non-periodic systems present characteris- 
tics of their own. We further discuss how most of our 
results can be understood within the framework of the 
collective coordinate technique || (see also the review 
H and references therein). Notwithstanding that ana- 
lytical insight, we have also found effects that cannot be 
interpreted in terms of such a particle-like behavior, and 
we have been able to associate those to the long-range 
order characteristics of aperiodic chains. 

The model we use as our working example is a damped, 
dc driven, discrete SG equation given by 
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(u n+ i - 2u n + u n -x) + V n smu n + au = F, (1) 



where dot means time derivative, a is the lattice spacing, 
and n runs over the lattice sites n = 1, . . . , N. The co- 
efficient in front of the on-site potential, V n , is directly 
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related to the physical properties of the application one 
is interested in: Thus, it has to do with local critical cur- 
rents in Josephson devices, or with the strength of hydro- 
gen bonds between complementary basis in DNA models. 
In the following, we will allow V n to take only on two val- 
ues, V a and Vb. Moreover, by an appropriate rescaling, 
it is possible to fix V a = 1, and so will be done hereafter. 
The spatial arrangement of the two kinds of values will 
be chosen to be either periodic, aperiodic, or random. 
As our aperiodic models, we pick two standard choices, 
namely the Fibonacci and the Thue-Morse chains. They 
are generated starting from two basic units A and B us- 
ing the following inflation rules: A — ► AB, B — > A for 
the Fibonacci chain and A — » AB, B — > BA for the 
Thue-Morse chain. In this way, finite and self-similar 
aperiodic chains are obtained by n successive applica- 
tions of these rules, with N = F n sites for the Fibonacci 
lattice and N — 2™ sites for the Thue-Morse lattice. Here 
F n = F n _i + F n _2 with F = F% = 1 are the Fibonacci 
numbers. The number of ^4-sites in the lattice is ~ tN 
in the Fibonacci case and N/2 in the Thue-Morse case, 
where r = hm n _ >00 (.F n _i/.F' n ) = (v5— 1)/2 is the inverse 
golden mean. Both chains have been used in very many 
contexts to model aperiodic ordering which, in spite of 
previous, more naive ideas, it is not something intermedi- 
ate between periodic and random systems (see, e.g., jlO| ] 
and references therein). 
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FIG. 1. Steady state velocity versus applied force for 
chains with different orderings and 14 = 10. Points cor- 
respond to numerical results, solid lines to the theoretical 



prediction in Eq. 



with V a vg corresponding to c = 0.5. 



(O), periodic chain with N = 4182; (□), Fibonacci chain 
with N = Fig = 4181; (A), Thue-Morse chain with 
TV = 2 12 = 4 096. Dashed lines joining symbols are a guide to 
the eye. The threshold for random chains is out of the range 
of the plot and its value is about F = 0.5. Inset: same but 
for Fibonacci chains with V a — 2, V a = 5, and V a = 10, from 
top to bottom. 

To characterize the dynamics of kink solitons on these 
systems, we have numerically simulated Eq. ([!]) by means 
of a fifth-order adaptive-stepsize Rungc-Kutta routine 
jllj which has been shown to be an accurate procedure 



(see and references therein). In the homogeneous 
case (V n = V a = 1), it has long been known Q that if a 
soliton initially at rest evolves according to Eq. ([[]) with 
F > 0, it eventually reaches a steady state in which it 
propagates along the chain with velocity 
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with the opposite sign if F < 0. This result was found 
by means of an standard collective coordinate calcula- 
tion, and the derivation can be found in 0. As another 
check of our simulations, we compared their outcome to 
this prediction and found an agreement better than 1% 
for all studied cases. We will also make use of this expres- 
sion, although in a modified form: If one assumes V n is 
constant and given by the average of V a and Vb weighted 
by their concentrations, say V n = V avg — cV a + (1 — c)Vb 
(here cN is the number of A-sites of the chain), and re- 
peats the same calculation in ||, the predicted final ve- 
locity is 
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i.e., the asymptotic velocity is predicted to be smaller 
(higher) than that of the homogeneous model when Vb > 

V a (V b < V a ). 

The results of our numerical simulations are collected 
in Fig. [j], where we present the value of v as a function 
of the applied force F for the different orderings consid- 
ered, always with Vb = 10. The plotted value of v was 
obtained by starting the simulation for each F with a 
kink at rest in the middle of the chain, as given by the 
exact continuum solution of the SG equation, and letting 
it evolve while monitoring its velocity until it reached a 
constant value. In all cases it was verified that the soliton 
shape remained almost unaltered, which is the necessary 
condition for the concept of velocity to make sense. The 
simulation parameters are a — 0.1, which is a prototyp- 
ical value and whose only effect is to fix the force scale, 
and a = 0.1, a value which yields a discrete chain but 
still close to the continuum to avoid side effects induced 
by pinning due to the effective Peierls-Nabarro poten- 
tial Q . In Fig. I] it is easily appreciated the different 
behavior of the different kinds of ordering considered. 
The periodic chain with alternating V a and Vb is very 
much accurately described by our theoretical prediction 
in Eq. (||) with average velocity V avg corresponding to a 
concentration of c = 0.5, whereas it is seen that neither 
Fibonacci (with c = r) nor Thue-Morse chains obey that 
equation. The first discrepancy arises as the existence 
of a threshold force, F c , below which solitons are pinned 
and do not move. A comparison of such F c for the ape- 
riodic chains with two random chains (not shown in the 
figure) with the same concentration of V a (r = 0.618 . . . 
for the Fibonacci case and 0.5 for the Thue-Morse case) 
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leads to the conclusion that the threshold is about three 
times higher for the fully disordered chain than for the 
corresponding aperiodic one. Above threshold, the ap- 
proximation ([}]) fails also to accurately predict the value 
of v, overestimating it appreciably for those non-periodic 
cases. However, the overall behavior of the random and 
the aperiodic chains is different, in the sense that in the 
random case, above threshold the v — F curve is non- 
monotonous and the final velocity depends strongly on 
the particular realization of disorder. As for the depen- 
dence of the results on the value of Vb, it can be seen 
in the inset of Fig. [j] that, while keeping Vb > V a , the 
higher V& the larger F c for Fibonacci lattices. On the 
other hand, when Vb = 0.1, i.e., much smaller than V a , 
we have found that, contrary to the naive expectations, 
there is still a threshold force for the Fibonacci case, 
about F c = 10~ 3 , which we have verified is not related to 
Peierls-Nabarro pinning because the periodic chain shows 
no pinning for such value of the force. This is a striking 
result that should be compared to the random chain one, 
where it is found that F c ~ 0.04 
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FIG. 2. Steady state velocity versus applied force for chains 
with different orderings. Points correspond to numerical re- 
sults, solid lines to the theoretical prediction in Eq. (^) with 
V m g corresponding to c = r. (□), 6th approximation (13 
sites); (O), 8th approximation (34 sites); (A), 9th approxi- 
mation (55 sites); (v)> 11th approximation (144 sites); (Q), 
full Fibonacci chain. All lenghts are close to 4 180 sites. A 
soliton at rest spans about 50 sites for comparison. Dashed 
lines joining symbols are a guide to the eye. Labels indicate 
the order of the approximation to the Fibonacci chain as well 
as the full aperiodic sequence. 

To gain further insight into the pinning phenomenon, 
we have compared the results for the full aperiodic chains 
with those of periodic chains with unit cell formed by a 
shorter approximation to the corresponding aperiodic se- 
quence. Thus, for the Fibonacci chain, we have studied 
periodic chains with unit cells V a VbV a (the 3 th Fibonacci 
approximation), V a VbV a V a Vb (the 4th approximation), 
and so forth up to the 11th approximation (144 sites in 
the unit cell). The results are shown in Fig. The most 
remarkable conclusion that can be drawn from this plot is 



that for periodic chains with unit cells smaller than the 
soliton size (around 50 sites when at rest), the behav- 
ior above threshold is independent of the unit cell, and 
moreover, it is very well described by the average veloc- 
ity introduced in Eq. (j^). This is very important, since it 
implies the existence of influences coming from the long- 
range order of the full aperiodic chain which do not arise 
in periodic approximations unless the size of the unit cell 
is much larger than the soliton width, i.e., unless the soli- 
ton is unable to distinguish the unit cell from the whole 
chain. Another remark in order here is that the threshold 
for the different periodic approximations depends on the 
initial position of the soliton in the chain, changing up 
to a factor two for different positions, although keeping 
below than that of the full Fibonacci chain. Results for 
the Thue-Morse model are basically the same, although 
in this case shorter lengths are needed for the simulation 
to get close to that of the full aperiodic chain, and the 
threshold dependence on the size is nonmonotonous. We 
tentatively associate this to the fact that in the Thuc- 
Morse sequence B-sites may appear in pairs, contrary to 
the Fibonacci sequence, and therefore it is to be expected 
that their effect will be stronger on the soliton. 

The existence of a threshold force is clearly the main 
failure of the collective coordinate theory for both the 
aperiodic and the random chains and, therefore, we un- 
dertook the task of finding a better analytical descrip- 
tion. To this end, we followed the same approach of the 
work by Kivshar and Salerno , where they introduced 
an effective potential to account for their results on DNA 
promoter dynamics. The basic idea is similar to that of 
the collective coordinate technique, but they improve it 
by including the spatial ordering of the chain. We skip 
the details here, as the interested reader may find them in 
Refs. and quote only the final result: The effective 

potential seen by a soliton, initially at rest at a lattice 
site no , is given by 

w , n > _ J2 m ( v avg + K)[sech 2 (z m ) - sech 2 (^ ) )] 
2X; m sech 2 (z m ) 

_ 2F 1 [(tan-V m )- tan'V")] 
2 Em sech 2 (z m ) 

with z m = aVavg(m - sn) and Zm) = aVUgijn - sn ), 
and the sums run over the whole lattice. Finally, to in- 
clude soliton width effects, the potential in Eq. (||) is 
averaged in the interval of the lattice spanned by the 
soliton. We have to stress that this approach only ap- 
plies to the early stages of the problem, when the kink is 
at rest or beginning to move at a very slow speed. This 
is so because in deriving Eq. (Q) the dissipative term is 
not included. Therefore, this approach should be use- 
ful in predicting the threshold force although it certainly 
does not apply to the dynamics above threshold. As can 
be seen from Fig. [| the agreement is very good for the 
Fibonacci chain, and the same can be said about the 
Thue-Morse chain and the random case (not shown; in 
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addition, the random case depends strongly on the real- 
ization considered, which was to be expected in view of 
the simulations). For comparison, notice that the poten- 
tial minimum is always absent in periodic lattices (see 
inset of Fig. [|), in agreement with our numerical simula- 
tions where pinning is not observed in those lattices. In 
computing the average value of W(n, no) we have used 
a soliton width of 40 sites, which is of course quite arbi- 
trary. We have checked that variations of ±10 sites are 
not crucial for our results, which remain semiquantita- 
tively correct. Indeed, due to this freedom in the election 
of the soliton width, we have not pursued a better agree- 
ment, because it would be difficult to justify the choice of 
that value aside from the fact that it fitted the numerical 
simulations. 
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FIG. 3. Effective potentials for a kink soliton initially at 
rest in the center of a Fibonacci chain for forces F = 0.05 
(solid line), F = 0.1 (dashed line), and F = 0.2 (dot-dashed 
line). Inset: Same, but for the periodic chain. 



In conclusion, we have studied soliton propagation 
along binary lattices with different orderings. We have 
found that whereas solitons can propagate for any force 
on periodic systems, non-periodic ones exhibit a thresh- 
old value, i.e., forces larger than a value F c are needed 
in order to start propagation. We have been able to an- 
alytically explain that, as well as to characterize the dif- 
ferences that, in turn, separate the phenomenology of 
random and aperiodic chains. We have also found that 
long-range order effects arise when solitons are propagat- 
ing along the chain, which show up in a decreasing of the 
steady state velocity with respect to the theoretical ex- 
pectations for the periodic chains. These results can be 
useful in the context of Josephson devices, as they can 
be the basis for the design of new devices with specific 
properties (the value F c corresponds to a critical current 
for the device to start conducing). This is also important 
as our results can be checked in specifically built JJA's. 
On the other hand, they can also be relevant to DNA 
promoter dynamics. Indeed, our results show that the 
long-range correlation present in DNA due to the infor- 
mation it encodes makes soliton propagation easier than 
if it were purely random, in fact allowing their traveling 



along the chain at lower velocities. This would not be 
possible if the structure of the molecule would be ran- 
dom in view of our results. 

It is clear that these results just opened the door to 
the problem of soliton propagation in aperiodic systems. 
As we have mentioned along the paper, there are a num- 
ber of unsolved questions, like an analytical explanation 
of the soliton velocity in aperiodic chains, or how peri- 
odic approximations converge to the full aperiodic sys- 
tem. Besides that, the study of other aperiodic models 
would be helpful in order to clarify the generic properties 
exhibited by this kind of models. Work along these lines 
is in progress. 

We thank E. Macia for illuminating conversations on 
aperiodic systems and J. A. Cuesta, F. Falo, L. M. Fiona, 
P. J. Martinez, and P. Quemerais for helpful discussions, 
and L. Garcia for the use of his computer facilities. 
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